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Abstract. 

Two samples of O- and B-type stars and Cepheids with Hipparcos data 
have been used to characterize galactic rotation and spiral arm kinematics 
in the solar neighbourhood. An extensive set of simulations has been per- 
formed in order to assess the capabilities of the method and its sensitivity 
to sample errors and biases in the adopted parameters. 



1. Introduction 

Young stars are good tracers of the galactic spiral structure in the solar 
neighbourhood, but, unfortunately, the observational data available up to 
now have not allowed to reach conclusive results on the parameters describ- 
ing such structure. The accurate astrometric data provided by Hipparcos 
(ESA 1997) allow us to re-examine the kinematic behaviour of these stars 
and, even more, to evaluate the possibility to undertake this study using 
more sophisticated kinematic models. 

In a previous paper (Torra et al. 2000) we characterized the structure 
and kinematics of the Gould Belt system, establishing its boundary to a 
distance of 600 pc from the Sun. Beyond this distance, the kinematic be- 
haviour of young stars, in particular 0-B stars and Cepheids, has to reflect 
the smooth galactic potential, that is the galactic rotation and the kine- 
matics owing to the spiral structure. The study of the stellar velocity field 
outside the Gould Belt is undertaken in the present paper, considering both 
the current observational material - its errors and the sample's constraints 
- and, through realistic simulations, the robustness of the resolution pro- 
cess. In particular, for Cepheid stars, the influence of the period-luminosity 
relation (PLR) on the derived kinematic parameters is also evaluated. 
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2. The working samples 

2.1. SAMPLE OF O- AND B-TYPE STARS 

The construction of this sample using Hipparcos data is fully described 
in Fernandez (1998) and Torra et al. (2000). Due to the uncertainty on 
trigonometric parallax at large distances, these were computed following 
Crawford's (1978) photometric calibration. Preliminary results applying the 
LM method (Luri et al. 1996) to Hipparcos data indicate this calibration 
is good enough for our purposes (Jordi et al. 2000). Considering only those 
stars with an heliocentric distance 0.6 < R < 2 kpc, the sample with 
distances and proper motions contains 444 stars, whereas the subsample 
with available radial velocities (Grenier 1997) contains 304 stars. 

2.2. SAMPLE OF CEPHEID STARS 

The initial sample contains all the Hipparcos classical Cepheids. Astromet- 

ric data were taken from this catalogue and radial velocities from Pont 
et al. (1994, 1997). Two PLR have been used for distance computation: 
Luri's (1999) relation, which corresponds to the short cosmic distance scale 
(J14 = -1.08 - 2.72 log P), and Feast k Catchpole (1997) relation, which 
gives a long cosmic scale (My = —1.41 — 2.81 log P). In both cases, peri- 
ods were obtained from Hipparcos and individual reddenings from Fernie 
et al.'s (1995) compilation (continuous updating). Up to R = A kpc, the 
sample with known distances and proper motions contains 203 stars and 
the subsample with known radial velocities contains 179 stars. 

3. Kinematic model 

Our kinematic model considers the systematic contributions of solar mo- 
tion, differential galactic rotation (up to second-order approximation) and 
spiral arm kinematics (modeled in the frame of the Lin's theory). In ad- 
dition to the solar motion, we derived the first- and second-order terms of 
the galactic rotation curve, the phase of the spiral structure at the Sun's 
position, the amplitudes in the perturbation of the velocity in the galac- 
tocentric and tangential directions and an adimensional parameter which 
takes into account the difference in the velocity dispersion between the 
solar-type stars and the considered stars. There are four free parameters: 
the number of spiral arms of the galaxy (m), their pitch angle (i), the galac- 
tocentric distance of the Sun (ro©) and the circular velocity at the Sun's 
position {@{zuq)). A couple of values of these parameters were considered: 
m = 2 or 4, i = -6° or -12°, wq = 7.1 or 8.5 kpc and 6(^0) = 184 
or 220 km s~^. The kinematic parameters were derived via weighted least 
squares fit, following an iterative scheme. Three cases were studied, solving 
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the equations considering only radial velocities, only proper motions and 
the combined resolution (see Fernandez et al. 2001). 

4. Test of robustness: simulations 

Numerical simulations were performed in order to quantitatively evaluate 
the biases in the kinematic model parameters induced by, among others, 
the observational errors, the incompleteness of the sample, the rejection of 
high residual stars and the correlations among the parameters. Both O- 
B star and Cepheid samples were simulated considering the same spatial 
distribution than the real sample and a Schwarzschild distribution for the 
residual velocities. Solar motion, galactic rotation and spiral arm kinematics 
were introduced following our kinematic model. We assumed gaussian errors 
in distance, radial velocity and proper motions. 

Detailed conclusions of these simulations can be found in Fernandez 
et al. (2001). As main results, we conclude that the combined solution 
minimizes the biases induced by correlations. The observational errors in 
proper motion, radial velocities and distance produce biases less than 0.5 
km in the solar motion components and in the galactic rotation and 
spiral structure parameters. Due to the spatial distribution of the available 
data, the best distance intervals were found to be 0.6 < R < 2 kpc for 0-B 
stars and R < 4 kpc for Cepheids. 

5. Results and discussion 

5.1. GALACTIC ROTATION CURVE 

For the A Oort constant, a discrepancy between the different solutions has 
been found: ^ 13.6±0.7 km kpc"^ for 0-B stars and A^^^- ^ 16.7± 
0.6 km s-i kpc-\ ^g^^ ^ 15.0 ± 0.6 km s-^ kpc-^ for Cepheids (Luri's 
1999 and Feast & Catchpole's 1997 PLRs, respectively). This discrepancy 
has been extensively found in the literature (e.g. Hanson 1987, Metzger et 
al. 1997, Feast & Whitelock 1997 and Mishurov & Zenina 1999 obtained 
values between 11.3 and 18.8 km s~^ kpc~^) and, as pointed out by Oiling 
& Merrifield (1998), it can be produced by the different spatial distribution 
of the samples. In agreement with Feast et al. (1998), we found a very 
small second-order term of the galactic rotation, compatible with a linear 
rotation curve in the solar neighbourhood. 

5.2. GALACTIC SPIRAL STRUCTURE 

Using both 0-B and Cepheid stars data, we found an angular rotation ve- 
locity for the spiral pattern of Qp 30-35 km s~^ kpc~^, incompatible with 
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the classical value 13.5 km kpc~ proposed by Lin and collaborators 
in the later 60s. Our value is in good agreement with the larger values re- 
cently published by Amaral & Lcpinc (1997) and Mishurov ct al. (1997, 
1999). Prom this value we can conclude that the Sun is outside - though 
near - the corotation circle. 

From the available data wc are not able to decide on a galaxy with 2 
or 4 spiral arms. Nowadays, this is an unsolved problem (see contradictory 
results in Amaral &: Lepine 1997, Mishurov &: Zenina 1999, Drimel 2000, 
Lepine et al. 2000). Nevertheless, good coherence is obtained to derive the 
phase of the spiral structure: from both samples, we obtained values in the 
rank [330°, 50°], indicating the Sun is near the potential minimum of an 
arm. The values of the velocity amplitudes due to the spiral perturbation 
were found to be small, about 2 km s"^. 
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